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ABSTRACT 

The recent detection of interstellar polarization in the solid CO feature near 
4.67/im shows that CO-mantled grains can be aligned in cold molecular clouds. 
These observations conflict with a theory of grain alignment which attributes 
the polarization in molecular clouds to the effects of cosmic rays: according to 
this theory, oblate spheroidal grains with H2O- and COvdominated ice mantles 
are spun up to suprathermal energies by molecular evaporation from cosmic ray 
impact sites but spin up does not occur for CO-mantled grains. Motivated by this 
conflict, we reexamine the effects of cosmic rays on the alignment of icy grains. 
We show that the systematic torques produced by cosmic rays are insufficient 
to cause suprathermal spin. In principle, the random torques due to cosmic 
rays can enhance the efficiency of Davis- Greenstein alignment by raising the 
grain rotational temperature. However, a significant enhancement would require 
cosmic ray fluxes 6-7 orders of magnitude larger than the flux in a typical cold 
cloud. 



Subject headings: magnetic fields - polarization - dust extinction 



Introduction 



The recent detection of polarization in the 4.67 /im feature of solid CO (Chrysostomou 
et al. 1996) sets severe constraints on theories of interstellar grain alignment. Solid CO 
mantles are expected to survive only in cold, dense clouds, where the temperatures of the 
gas (T g ) and grain solid material (T^) must be nearly equal due to collisional coupling. For 
example, the Davis- Greenstein alignment mechanism is a dissipative process which is driven 
by the difference between T g and Td (Jones & Spitzer 1967). Recent calculations (DeGraff, 
Roberge & Flaherty 1997) show that the largest polarizations observed toward molecular 
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clouds are inconsistent with Davis-Greenstein alignment unless T g /T^ ^ 10 and that this 
conclusion holds even if the grains are superparamagnetic. [| Since a temperature disparity 
of this magnitude is inconsistent with observations, it seems certain that the grains in cold 
clouds — including CO-mantled grains — are aligned by some "nonthermal" process. 

The first nonthermal model of grain alignment was developed by Purcell (1975, 1979, 
hereafter P79), who pointed out that the temperatures of the gas and dust do not limit the 
efficiency of alignment if the grain rotational energies are sufficiently large. P79 pointed out 
that spin up to suprathermal energies will occur whenever the torque on a grain, referred to 
axes fixed in the grain material, has a nonvanishing time average over times ^ the timescale 
for frictional coupling to the gas. P79 also identified three mechanisms that can produce 
such a "pinwheel torque": (i) the formation of H 2 at catalytic sites on the grain surface, 
which requires atomic hydrogen; (ii) gas-grain collisions on a surface with spatial variations 
in the "accommodation coefficient," which requires T g ^ Ta; and (iii) photoelectric emission 
from a surface with variations in the photoelectric yield, which requires UV photons. All of 
these processes are suppressed in dark clouds, where the atomic hydrogen concentration and 
UV flux are negligible and T g « T^. If the starlight flux in the blue part of the spectrum^ is 
negligible, then we may also disregard the radiative torques arising from differential scattering 
(Dolginov & Mytrophanov 1975, Lazarian 1995a, Draine & Weingartner 1996a,b). 

In view of these circumstances, it is natural to inquire whether some other process might 
cause suprathermal rotation. A natural possibility is related to the evaporation of molecules 
adsorbed on grain mantles: if a grain is not heated uniformly, then hotter places on the 
mantle will desorb molecules at higher rates compared to colder places. These "hotspots" 
can act in a manner that is similar to the catalytic sites of H 2 formation in Purcell's model 
for suprathermal spin up. The required temperature nonuniformity might be caused by 
nonuniformities in the absorption of light. Consider, for example, a small soot inclusion on 
the grain surface. Such an inclusion would absorb light and have a temperature greater than 
the mean temperature of the grain. In other circumstances, a darker inclusion would radiate 
more efficiently and therefore become colder than the rest of the grain. Since variations in the 
grain absorption (emissivity) are likely to persist for much more than a gas damping time, 
the resulting torques should be long lived. If the grain dynamics are dominated by these 



1 However, helical grains may be aligned even if T g w (see Lazarian 1995). The process is discussed at 
length in Lazarian, Goodman & Myers (1996). 

2 For "standard" grains of size ~ 10 _5 cm, only short-wavelength radiation causes an appreciable torque 
(Draine & Weingartner 1996a). However, if the aligned grains in molecular clouds are substantially larger, 
then longer wavelengths will affect the alignment. 
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torques, rather than the torques due to gas damping^ the alignment could be substantial. 
We plan to address this interesting possibility in a subsequent paper. 

In this paper, we examine the feasibility of a related mechanism suggested by Sorrell 
(1995a,b). Sorrell pointed out that a cosmic ray impact heats a grain locally, and that 
molecular evaporation from the resulting hot spot will cause rotational acceleration via the 
rocket effect. According to Sorrell's analysis, the rocket effect on an individual grain produces 
a nonzero time-average torque with a correlation time that is limited by changes in the mantle 
surface. In this view, the rocket effect spins up the grain to a suprathermal kinetic energy 
and the alignment occurs via Purcell's mechanism (P79). However, the scenario considered 
by Sorrell (1995a) conflicts with the observations of Chrysostomou et al. (1996): Sorrell's 
model predicts that grains with H 2 0- and C0 2 -dominated mantles will spin up to energies 
~ 100/cTg but that spin up of CO-mantled grains does not occur. In an effort to understand 
this conflict, we reexamine the effects of cosmic rays on the alignment of ice-mantled grains. 
In §2, we model the evaporation of molecules from a cosmic ray hotspot. In §3, we consider 
the conditions under which evaporating molecules may produce a pinwheel torque and hence 
suprathermal rotation. In §4, we describe the possible effects of the random torques due to 
evaporation and the resulting enhancement of Davis- Greenstein alignment. We summarise 
our results in 55. 



2. Evaporation from a cosmic ray hot spot 

The evaporation of molecules caused by cosmic ray heating has been discussed elsewhere 
(Watson & Salpeter 1972; de Jong & Kamijo 1973; Aannestad & Kenyon 1979; Leger, Jura 
& Omont 1985, henceforth LJ085). Here we merely apply the results of LJ085 to model 
the evaporation of H 2 and CO molecules from H20-ice mantles. A cosmic ray nucleus with 
charge Ze and energy E loses energy in solid H 2 at a rate 



Q(Z,E) 



9.40xlO- 4 (|) Z (gS,)~°- 75 erg cm" 1 



0.02 < < 0.2 

GcV 



1.48 xlO- 3 (|) 2 [l + O.l (^) 



-1.5 



(2-1) 



erg cm" 1 , 0.2 < ^ < 10 



3 For this purpose, the randomization of angular momentum during a spin-up interval should be small 
(Lazarian 1995b). 
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(LJ085).[] The energy is deposited initially (on a timescale ~ 10~ n s) in a cylinder of radius 
r ~ 50 A around the cosmic ray track and spreads laterally thereafter by thermal diffusion. 
Over the temperature range of interest here, the thermal diffusivity of H 2 ice depends only 
weakly on temperature (Zeller & Pohl 1971). Consequently, the heating of the mantle is well 
described by the heat diffusion equation with constant thermal diffusivity. The solution for 
the energy density, U, is 



U(r,t) 



Q 



Ana (t + 1 ) 



exp 



4a (t + 1 ) 



(2-2) 



where r is distance perpendicular to the cosmic ray track, t is time, a is the thermal diffu- 
sivity, and 

Here we have adopted somewhat arbitrary initial conditions, such that U(r, 0) is a Gaussian 
function of r with a fraction / of the total energy contained in r < r at t = 0. In the 
following discussion, we will assume that r = 50 A and set / = 0.5. The temperature, T, 
can be found from the relation 

r T(r,t) 

P C v (T>)dT>, (2-4) 



U(r,t) 



where To is the temperature before the impact and pCy is the volume specific heat. In 
the calculations discussed below, we set T = 15 K and adopt the values of a and pCy 
recommended by LJ085 for H 2 ice. 



We assume that molecules evaporate from a hotspot at the classical rate, 

R vap (r,t) = vo exp [-AH s /kT(r,t)) , 



(2-5) 



where uq and AH S are respectively the lattice vibration frequency and binding energy of a 
molecule adsorbed on the mantle surface. The probability that a molecule at radius r has 
evaporated before time t is therefore 

ft 



Pnp(r,t) = 1 - exp - J 



^vapM') dt ' 



(2-6) 



The mean value of the total number of molecules evaporated by a single cosmic ray with 
charge Z and energy E is 

N vav (Z,E) = A ea (Z,E)6 s , (2-7) 



4 The energy loss rates in eq. (2-1) include the corrections derived by LJ085 to account for oblique cosmic 
ray impacts and partial escape from the grain of electrons ejected by the cosmic ray. We have taken the 
density of H2O ice to be 1 gem -3 . 
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where 8 S ( cm 2 ) is the surface coverage before the impact and 

A eS (Z,E)= rfr27rrP vap (r,t max ) (2-8) 
J o 

is the effective area "cleaned off" by the cosmic ray. We assume that no evaporation occurs 
for times t > t max , where t max (r) is the time when the mantle at radius r has cooled to 
temperature T ; because virtually no evaporation occurs for T « T , the arbitrariness of this 
definition has virtually no effect on our numerical results. The value of A e g depends on Z 
and E only through the energy loss rate, Q. The dependence of A e s on Q is illustrated in 
Figure 1 for the evaporation of H2 (^0 — 7.5 x 10 12 , AH S / k = 555 K; Sandford & Allamandola 
1993) and CO (z/ = 7.0 x 10 14 , AHjk = 1030 K; LJ085) molecules from H 2 ice. 

The temperature of the vapor at a hot spot 

rp /y p\ — /o°° dr2nr / °° dt [1 - P vap (r,t)] R vap (r,t) T(r,t) 
Vapl ' } ~ I0 00 dr 2vrr / °° dt [1 - P vap (r, t)] P vap (r, t) ^' J> 

is the mean (for all evaporating molecules) temperature of the surface at the instant of 
evaporation. The numerical evaluation of expression Q2-9|) for CO and H 2 evaporations 
shows that, for both H2 and CO, T vap increases from w 50 K to ^ 200 K as Q increases from 
10~ 4 ergcm _1 to 10 _2 ergcm _1 . 



3. Constraints on the pinwheel torque due to cosmic rays 

The evaporation of a molecule from the mantle surface produces an impulsive change 
in the grain angular momentum, J, due to the rocket effect. The rotational dynamics of a 
grain subject to many such impulses are determined by the mean torque, 

and diffusion tensor, 

/AJAJ\ . . 

where AJ is the cumulative change in J caused by impacts during a time interval At and 
the angle brackets denote time averages. The mean torque due to the rocket effect has 
components 

4 Ck) = ^dE r dE ME) iV vap (Z, E) JT k , (3-3) 

where Sd is the grain surface area, 4>z(E) dE is the omnidirectional flux of cosmic rays with 
charge Ze and energies between E and E + dE, and 5Jk is the mean angular impulse due 
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to a single evaporating molecule. In deriving expression (|3-3|) , we have assumed that the 
flux of cosmic rays is isotropic in the grain frame and that each cosmic ray impact creates 2 
hotspots. The analogous expression for the diffusion tensor is 



B ki 



1 roo 

-S d ^2 dE<j> z (E)N yap (Z,E)5J k 8J l . (3-4) 



(rck) _ 

~~ 9 ^ d I 

Z 



The quantities SJk and 5JkSJi are determined by averaging the angular impulse due to a 
single evaporation over the position of the evaporation site and momentum of the evaporating 
molecule. Consequently, 5 Jk depends in general on the grain shape. In Appendix A, we show 
that 5Jk = for an oblate spheroid, the shape adopted by Sorrell (1995a,b). It follows that 

(rck) 

A k =0, that is, there is no pinwheel torque due to cosmic ray impacts on an oblate 
spheroid. This statement is true for any surface of revolution. 

Of course, real interstellar grains are not surfaces of revolution and it is possible to show 
that pinwhell torques do not vanish for less symmetric shapes (e.g. they do not vanish for 
square prisms). Nevertheless, we can place an upper limit on the magnitude of the pinwheel 
torque for any shape by assuming, unrealistically, that all of the angular impulses produced 
by evaporating molecules lie along the same direction in the grain frame. Since the cosmic 
ray hits on a real grain are uniformly distributed over its surface, this assumption is obviously 
extremely optimistic. 

Let 5 J = b (mvapfcTvap) 1 / 2 be the mean magnitude of an individual impulse, where b is 
some characteristic linear dimension of our hypothetical grain. In our optimistic scenario, 
the mean torque due to the rocket effect would have magnitude 

1 POO 

A^ = -S d 6 s ^ / dE<f) Z (E)A CTl (Z,E)5J(E). (3-5) 

Now suppose that all of the angular impulses lie along the kth principal axis of inertia. If we 
assume that the rotational friction is provided by gas damping, then the grain would spin 
up to kinetic energy E TOt , such that 



E, 



rot 



A (rck) t 



gas,fc 



\kT g hkTg 



(3-6) 



where Ik and t gas ,A: are respectively the rotational inertia and gas damping time for rotation 
about axis k. 



We have evaluated expression (|3-6|) using the functional form of 4>z given by LJ085, 
cosmic ray abundances given in Simpson (1983) for the local ISM, and typical values for 
the grain properties (b = 10~ 5 cm, T d = 15 K, p d = 3 gem -3 ) and physical conditions 
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(n g = 10 4 cm -3 , T g = 15 K) in the gas. We find the same value, E rot ~ 10 -5 kT g , whether 
we assume that the evaporating molecule is H 2 or CO. Note that this is the average excess 
energy for an ensemble of grains. For the cosmic ray fluxes adopted here, most of the grains 
in the ensemble actually experience no cosmic ray hits during the interval t ga , S} k, so the 
average energy is somewhat misleading. However, we estimate that the largest energies in 
the ensemble are only ~ 10 _3 kT g , so this does not alter our conclusion. 

Evidently, cosmic rays cannot spin up the grains under the most optimistic scenario that 
one can conceive. It is easy to understand this result if one notes that only heavy cosmic 
rays have energy loss rates sufficiently large to produce a significant number of evaporations. 
For example, the cosmic ray spectrum calculated by LJ085 peaks at E « 0.3 GeV for both 
protons and iron nuclei. At E = 0.3 GeV, the energy loss rates for protons and iron nuclei 
are Q p = 4 x 10~ 6 erg cm -1 and Qp e = 2 x 10 -3 erg cm -1 , respectively. According to Figure 1, 
evaporations caused by protons impacts would be completely negligible at this energy. 



4. Rotational excitation by random cosmic ray torques 



The stochastic torque produced by the rocket effect can enhance Davis-Greenstein align- 
ment by increasing the grain rotational temperature (Salpeter & Wickramsinghe 1969; Pur- 
cell & Spitzer 1971). The random angular impulses produced by evaporating molecules cause 
the angular momentum to change in random walk fashion so that, in the absence of other 



processes, the kth component of J would increase without limit as Jk oc 



1/2 



In re- 



ality, the rotational friction produced by gas damping and other dissipative processes limits 

r (rck) 1 l/^ 

the growth to Jk ~ B\ k tdamp , where tdamp is the relevant damping time. For exam- 
ple, suppose that the only interactions of the grains with their environment are provided by 
gas damping plus the rocket effect. Then one can show (see RDGF93, eq. [3.18]) that the 
distribution of Jk is Maxwellian with an effective temperature 



T, 



eff 



n kk + D 



(rck) 
kk 



^gas,k 



2I,k 



(4-1) 



where Bf k is the diffusion tensor for gas damping. If B kk ^ is sufficiently large, then the 
rotational excitation provided by the cosmic rays can permit Davis-Greenstein alignment 
even in clouds with Td = T„. 



It follows from equation ( [4-l| ) that the stochastic torque produced by the rocket effect 
increases the grain rotational temperature by a factor 1 + 0, where 



= 



B 



(rck) 



kk 



B 



(g) 

kk 



(4-2) 
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After combining expressions ( |2-7|) , (3-4), and (A2), we find that the nonzero components of 
the diffusion tensor for the rocket effect are 



6 z 



Air 



where 



W z = dE d> z (E) A eS (Z, E) T^Z, E) 



(4-3) 



(4-4) 



and the Tk are weak functions of the grain eccentricity with 3/8 < < 1 (see Appendix A). 
The diffusion tensor for gas damping is also diagonal with components 



( g ) ™2,4 ,3 ( -. , Td 

B kk - 1 k n g m g b v th I 1 + — 



(4-5) 



v t h = \/2kTg/m g is the gas thermal speed. It follows that 



(RDGF93), where n g is the number density of the gas, m g is the mass of a gas particle, and 



(4-6) 



= v^(l + T d /T g 



m 



vap 



EzW z 



Notice that is independent of the surface eccentricity. 

In order to estimate 0, we used the results of §2 to evaluate expression ( |4-4p by numerical 
integration. Due to the steep dependence of A efi on Q (Fig. 1), the sum is dominated by iron, 
the most abundant nucleus with Z> 1. For example, consider the relative contributions of 
iron and cosmic ray protons. LJ085 found that iron and protons have approximately the 
same energy spectrum; according to their models, <pz = AttAzI, where Az is the abundance 
of cosmic rays with charge Ze relative to the abundance of protons and 



1(E) 



r 2i (civ) 

1.5 

0-3 (c&) 
0-3 (M 



-i 



0.02 < < 0.07 
0.07 < ^ < 0.2 
2 < < 1 

U ' Z ^ GeV ^ 1 



(4-7) 



1 < 



GeV 



(cm -2 s" 1 GeV -1 ) is the mean intensity of cosmic ray protons. Taking A26 = 3 x 10 -5 , the 
cosmic ray abundance of iron in the local interstellar medium (Simpson 1983 and references 
therein), we find that W2& = 5 x 10~ 13 K s _1 and W\ = 4 x 10~ 14 K s _1 for the evaporation of 
H2 from H2O ice. The analogous calculations for CO evaporation yield W26 = 3 x 10~ 13 K s _1 
and W\ = 6 x 10~ 17 K s -1 . We will assume henceforth that all of the evaporation is caused 
by iron cosmic rays and set J2 Z Wz = 4 x 10~ 14 K s _1 for H 2 evaporation and J2 Z Wz = 
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5 x 1CT 13 K s -1 for H 2 evaporation and J2zW z = 3 x 1CT 13 K s -1 for CO evaporation.^ To 
estimate O for typical cold cloud conditions, we will assume that the gas is composed of 
H 2 , with n g = 10 4 cm~ 3 , T g = T d = 15 K, and 9 S = 10 15 cm~ 2 . With these assumptions, we 
find that © ~ 10~ 7 if the rocket effect is caused by the evaporation of H 2 and ~ 10~ 6 if 
the evaporating molecules are CO. We conclude that the rotational excitation provided by 
the rocket effect is insignificant unless the flux of iron nuclei were increased by 6-7 orders of 
magnitude. However, an increase in <p of this magnitude can be ruled out in general: with 
the flux adopted in our estimates, iron nuclei would already contribute about 10% of the 
total cosmic ray ionization rate, (, in a typical cloud with ( ~ 10~ 17 s _1 (LJ085). 

In principle, the desorption produced by low-energy electrons (Johnson 1990) could also 
contribute to the rotational excitation. However, the required electron flux is ruled out 
unambiguously by data on the interstellar ionization. 

Purcell & Spitzer (1971) computed the rotational excitation which is associated with 
energy loss by cosmic rays in the grain solid material. Comparing our results with those 
of Purcell & Spitzer (1971), we conclude that energy loss is more important for enhancing 
the rotational temperature of grains than the desorption of molecules from hotspots. Nev- 
ertheless, even the process studied by Purcell & Spitzer has only a marginal effect on grain 
alignment. 



5. Summary 

Our study has shown that the evaporation of molecules from cosmic ray hotspots pro- 
duces no pinwheel torque on an oblate spheroidal grain, regardless of its mantle composi- 
tion. Although it is possible, in principle, that nonzero pinwheel torques exist for some grain 
shapes, we have shown that the resulting increase in the grain rotational energy is negligible 
for realistic fluxes of heavy cosmic rays. Our findings resolve the apparent conflict between 
the recent observations of Chrysostomou et al. (1996) and the grain alignment model of 
Sorrell (1995a), inasmuch as the alignment process postulated in Sorrell's model does not 
occur. We have also shown that the random torques caused by cosmic ray hotspots produce 
rotational excitation that can enhance the efficiency of Davis-Greenstein alignment, even in 



5 To estimate the constributions of other nuclei, we evaluated J^z ^z by setting <f>z = 4<kAzI, assuming 
that I is given by expression (4-7), and adopting cosmic ray abundances appropriate for the local ISM 
(Simpson 1983). In this approximation, adding the effects of nuclei with Z ^ 26 increases J2 Z Wz by less 
than a factor of two for H2 and CO evaporation. An error of this magnitude in Wz has no effect on our 
conclusions. 
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clouds where = T g . However, the enhancement is completely insignificant for reasonable 
estimates of the cosmic ray flux and less than the enhancement associated with cosmic ray 
energy losses in the grain material (Purcell & Spitzer 1971). 

We are grateful to Bruce Draine for valuable comments. A.L. acknowledges the support 
of NASA grant NAG 5 2858 and W.G.R. acknowledges the support of NASA grant NAGW- 
3001. 



A. The quantities 5Jk and 5Jk5Ji for an oblate spheroid 

In order to compare our results with those of Sorrell (1995a), we assume that the mantle 
surface is an oblate spheroid. If the cosmic ray flux is isotropic in the grain frame, then the 
evaporation sites are uniformly distributed over the grain surface. The momentum distri- 
bution of the evaporating molecules is determined for thermal evaporation by the principle 
of detailed balancing (e.g., Roberge et al. 1993, hereafter RDGF93). Then the required 
coefficients are[] 

5J k = 0, k = x,y,z, (Al) 

and 

SJ k SJi = m vap /cT vap 6 2 e fc 5 H k = x,y,z (A2) 

(RDGF93, cf. eqs. [B8] and [BIO]), where h is the grain radius and the components are 
relative to a Cartesian basis with z parallel to the grain symmetry axis. The quantities 

£fc(e) = — ; ; — ^ ~L — ttt> k — x,y,z, (A3) 

V ; 3[1 + (1 -e 2 )g(e)Y ' y ' ' K ' 

are weak functions of the mantle shape with l/2<efe<2/3 and the geometrical factors 

Tz = TE i 3 + 4(1 " e2)9{e) " e_2 1 1 " (1 " e2)2g{e) \ } ' (A4) 



In eqs. [Al]— [A2], we have neglected terms of order jj^—, where ft is the grain angular velocity and 



iVap is a typical thermal velocity for an evaporating molecule. Here we have anticipated the result (§3) that 
the rotational energies of the grains are nearly thermal, so that flb/v v&p ~ \/m vap /Ma, where m va p is the 
mass of an evaporating molecule and Ma is the grain mass. Note that, although SJk — to zeroth order 
in rife/w vap , the first-order term we have neglected is always negative, so that evaporation from cosmic ray 
hotspots would actually spin down an oblate spheroid. We will neglect the rotational damping caused by 
cosmic rays, which is much smaller than gas damping for reasonable estimates of the cosmic ray flux. 
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and 



r, = r, = -{7-e 2 + (i 

where 



e 2 y 9 (e) + (1 - 2e 2 ) 



1 + e 



f \ 1 l /1 + e 
9(e) = — In 



2e VI -e 



l-(l-e 2 )^(e)]]}, (A5) 

(A6) 



are weak functions of e with 3/8 < F k < 1. 
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FIGURE CAPTION 

Fig. 1 — The mean number of molecules evaporated by a cosmic ray with energy loss rate 
Q equals the number of molecules initially adsorbed in the effective area A ef[ (Z, E) (see eq. 
|2-8|| ) . Results are shown for the evaporation of H 2 (solid curves) and CO (dashed curves). 
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